Slow interaction ramps in trapped many-particle systems: universal deviations from adiabaticity 
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For harmonic-trapped atomic systems, we report universal non-adiabatic features in the response to interac- 
tion ramps. We provide results for several different systems in one, two, and three dimensions: bosonic and 
fermionic Hubbard models realized through optical lattices, a Bose-Einstein condensate, a fermionic superfluid 
and a fermi liquid. The deviation from adiabaticity is characterized through the heating or excitation energy 
produced during the ramp. We find that the dependence of the heat on the ramp time is sensitive to the ramp 
protocol but has aspects unexpectedly common to all systems considered. We explain these common features in 
terms of universal dynamics of the system size or cloud radius. 
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Introduction — Adiabaticity is an essential and ubiquitous 
concept in quantum dynamics. In the current era of many 
novel non-equilibrium experimental possibilities, deviations 
from adiabaticity in slow parameter changes have attracted a 
lot of attention ifl Uloll . The question of non-adiabaticity is 
of fundamental interest, but also has practical implications. 
Many experimental protocols involve adiabatically changing 
a parameter in order to reach a desired quantum state. Since 
non-adiabatic heating can rarely be completely avoided, it is 
essential to understand deviations from adiabaticity in slow 
ramps. While the effect of quantum critical points in the 
ramp path has been considered in much detail ||il|2|], settings 
for non-equilibrium experiments in isolated systems tend to 
be mesoscopic rather than macroscopic, without true quan- 
tum critical points. Understanding non-adiabatic ramps in 
finite quantum systems is therefore vital. Also, since cold 
atoms dominate experimental non-equilibrium studies, a har- 
monically trapped many-particle system is the most important 
paradigm today for studying quenches and ramps. A few stud- 
ies of ramps in finite and trap ped systems have appeared in the 
very recent literature ioi llOll . indicating an emerging recogni- 
tion of the importance of the adiabaticity issue in finite sys- 
tems. In addition, in the past couple of years, reports of ex- 
perimental investigations of finite-rate interaction ramps have 
started to a ppe ar, both in the continuum jl III and in optical- 
lattices lll2l4l4ll . with harmonically trapped atoms. 




In this work, we consider non-adiabatic ramps in sev- 
eral distinct interacting many-particle systems, confined in 
isotropic harmonic traps. In each case, we consider ramps 
of the interaction from an initial value Ui to a final value Ui, 
occurring in time scale t. We focus on large but finite t 
(near-adiabatic ramps). We study deviations from adiabaticity 
through the heating Q, which is the final energy at time t ^ t 
minus the ground state energy of the final Hamiltonian. This 
quantity is also called the residual energy or excess excitation 
energy IslSlstl, and may be thought of as the "friction" due to 
imperfect adiabaticity lisll . The asymptotic form of Q{t) is 
a quantitative characterization of minimal corrections to adi- 
abaticity. Ref. |fl4ll reports measurements of excess energies 
after ramps, which may be regarded as the first experimental 
approach to this quantity. 

We find that the asymptotics of Q{t) is common to all the 
compressible systems that we considered: we have thus un- 
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FIG. 1: Ramp shapes. We illustrate Ui > Ui; it should be straight- 
forward to visualize U{ < Ui cases. 



covered universal non-adiabatic features of ramp dynamics in 
trapped systems. The Q{t) function has overall power-law 
decay, Q ^ t~'^ , with the exponent v depending on the shape 
of the ramp. For certain ramp shapes, Q{t) has oscillations 
superposed on top of the power-law decay. We present re- 
sults for a range of systems, interactions, and dimensional- 
ities, which make clear that the universal features are inde- 
pendent of system details and generally do not depend on the 
initial and final values of the interaction, as long as the trapped 
system remains in the same phase. 

Since the effects are universal over a wide range of trapped 
systems, they should be due to some type of dynamics that 
is prevalent in many harmonically trapped systems. We show 
that the relevant dynamics is the size oscillation or breathing- 
mode oscillation of the trapped cloud. Almost all trapped sys- 
tems have "soft" breathing modes due to vanishing density at 
the edge, irrespective of the nature of intrinsic modes of the 
system. Our results show that, in the slow ramp limit, the 
breathing modes due to trapping dominate the near-adiabatic 
response of many systems. To show that the asymptotic be- 
haviors of Q{t) are due to size dynamics, we will use a vari- 
ational description for one of the systems (the Bose conden- 
sate), treating the extent (radius) of the many-particle cloud 
as a time-dependent variational parameter. Such a "radius dy- 
namics" description will be shown to reproduce the universal 
Q{t) behaviors in thorough detail. An equivalent formulation 
is not easy to set up for all the systems; nevertheless, the com- 
monality of the Q{t) features, and the success of the radius 
description for at least one case, is convincing argument that 
the same dynamics type is the relevant feature in each case. 

Trapped atoms are by far the most promising setup for ex- 
perimentally exploring isolated-system dynamics in general 
and non-adiabaticiy issues in particular. A universal excita- 
tion mechanism that is dominant in generic trapped systems 
is thus an important baseline perspective for understanding 
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the many further non-adiabatic ramp experiments exptected 
in the near future. Contemporary theoretical treatments of 
non-adiabaticity in many-particle systems almost invariably 
appeal to cold-atom experiments for motivation. Yet, our re- 
sults show that in a real trapped-atom experiment radius dy- 
namics dominates over the intrinsic heating mechanisms that 
may be important in individual uniform systems. 

Distinct trapped systems — We will present results for 
fermionic and bosonic systems, with and without optical 
lattices. The lattice systems are two-component fermions, 
and single-component bosons, described respectively by the 
fermionic and bosonic Hubbard models: 

<ij>(j j 
<ij> j 

We use the (inverse) hopping strength as energy (time) units; 
the notations are otherwise standard. In one dimension (ID), 
the trap terms are H*^" — \kt^^pc}-^Cja for fermions 

and \kt^^pb^jbj for bosons. We will consider both spin- 
balanced {N^ = N]) and polarized (N^^N^) fermionic sys- 
tems. We treat time evolution of the lattice systems through 
exact numerical evolution of the full quantum wavefunc- 
tion, for few-particle configurations in ID. The Bose-Hubbard 
model is also treated via the Gutzwiller approximation, for 
iV > 10 bosons and/or for 2D. 

The continuum systems (without optical lattice) are the 
Bose-Einstein condensate (BEC) and the interacting two- 
component Fermi gas. We use trap units for these cases, ex- 
pressing lengths (energies) in units of trap oscillator length 
(trapping frequency). The BEC is treated via the Gross- 
Pi taevskii (GP) description lfl6l[l7ll . The dynamics is given by 
the time-dependent GP equation, = — ^ V' + \'r'^'4' + 
U {t)\ip\'^tlj, and the GP energy functional is 

EH'] = I + \u{im^ + ir^iv^p] . 

J r 

Here r is the radial position; = Jd^r is the spatial integral 
for dimensionality D; U is the effective interaction strength 
whose relation to the physical interaction is also D-dependent 
(c.f. Ref. in for ID). We normalize |i/f = 1, so that U 
contains a factor of the boson number N; thus U can be large 
even within the mean-field regime. We will present GP results 
general to all D. 

In addition to full solutions of the GP equation, we also use 
a variational description which uses the radius or size of the 
condensate, a{t), as the time-dependent parameter. Using a 
gaussian ansatz for the cloud shape, the equation of motion 
for a is found to be 

aa" + - a-^ - {V2^(7)-"U{t) =0, (1) 

the primes denoting time derivatives. The energy is 

E[a] = \D + ^2 ^ (^a'f] + i(V2^a)"^C/(t) . (2) 



This radius description is suitable for describing breathing- 
mode oscillations. For constant U, for small-amplitude oscil- 
latory solutions (7 = i?o + psin(f2t) (amplitude p), Eq. (|2|i 
shows that the excitation energy scales as ^p-^. We could also 
use a Thomas-Fermi instead of Gaussian profile; the results 
are very similar and do not affect our arguments. 

Finally, the continuum fermionic system (3D) will be 
treated in the hydrodynamic approximation. There are sev- 
eral similar formulations; we choose the so-called "time- 
dependent DFT" where the fermi gas is described by 
a nonlinear Schrodinger equation, = [— ^V^ + ^r^ + 
li{n, U)]^). Here n ~ is the sum of the (equal) densi- 
ties of the two components. We use the Hartree expression: 
p{n,U) = ^{Sn^n)"^/^ + \Un. Both a paired superfluid 
{U < 0) and a Fermi liquid {U > 0) are described by the 
same formalism. Analogous to the BEC case, we use a vari- 
ational description to formulate the dynamics in terms of the 
cloud radius. 

Ramp shapes — We analyze ramps of the form 

u{t) ^ Ui + e{t) {Uf - Ui) s{t/T) . 

The ramp function s{x) starts at .s(0) = and ends at 
s(oo) = 1. The ramps take place over time scale r but we 
do not require them to end at t = r. (Contrast, e.g., Ref. iQl ) 
We choose a collection of ramps which allows us to compare 
the presence/absence of kinks and various exponents. Specif- 
ically, we consider the following forms for s{x): 

[lA] X 61(1 - x) + e{x - 1) [IB] 1 - e"^ 
[2A] x^ e{l - x) + e{x - I) [2B] l-e"^' 
[3A] x^ e{l ~ x) + e{x - I) [3B] l-e~^' 

Each [A], [B] pair has the same initial behavior, s{x)^x'^, but 
the [B] versions have no endpoint kinks (Fig.[T]i. 

Universal features of ramp response — In Fig.|2]we present 
the behavior of the heat function Q{t), normalized against 
its instantaneous-quench value Q{t — 0). The small-r be- 
havior and the exact magnitude of Q{t) are system- and 
approximation-dependent; the universal features we present 
concern only the asymptotics of Q(r)/Q(0). 

In the left panels, the Q{t) behaviors are compared for the 
ramp shapes s{x)^x°' with discontinuous derivatives at end- 
points, [1 A], [2A], [3A]. Each curve has an overall power-law 
decay with the same decay exponent, Q{t) ^ t^^. This sug- 
gests that the residual energy for such ramps is primarily set 
by the endpoint kink. Superposed on the power-law decay are 
oscillations, which are often but not always smaller for larger 
a (contrast panels (c) with others). The oscillation strength 
decays faster for larger a, as r~^". 

The right panels concern smoothed ramps [IB], [2B], [3B], 
which lead to non-oscillating decay of Q(r). The decay expo- 
nent is seen to depend on the power a of s{x) ^ x", namely, 

Q(t) 

The dimensionality does not affect the decay exponents or 
general behavior. As an example, results are shown for the 
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[A] - type rampshapes 



[B] - type rampshapes 
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FIG. 2: The normalized excitation energy or heat, Q(r)/Q(0), for 
various trapped systems. (Ui,Uf) values and trap constants fctr are 
chosen over a wide range to emphasize the parameter-independence 
of the effects. Response to [A] ramps (left panels) shows overall 
Q{t) ~ T~'^ decay, with superposed oscillations that decay faster 
for larger-a ramps. For [B] ramps (right panels) the asymptotic be- 
havior is pure r~^° decay. Insets to (c): density profiles for U = f/i,f 
ground states. For (f), U = Ui^f ground-state profiles are very simi- 
lar, only Ui shown. 



BEC in ID, 2D, and 3D, in the (d) panels and insets. We have 
also found that the same exponents and oscillation features 
also appear in additional cases not shown, e.g., a continuum 
Fermi superfluid (3D, U < 0), the Bose-Hubbard model in 
higher D (treated via the Gutzwiller approximation), etc. 

The (c) and (f) panels involve systems which cannot be de- 
scribed as having single-radius profiles. The two spin com- 
ponents have different extents in the spin-imbalanced Hub- 



bard model (c); the Bose-Hubbard situation (f) has superfluid 
wings around a Mott core. The Q{t) behaviors are more rich 
for these systems; however, the features discussed above re- 
markably also persist in these more complex cases. 

Radius dynamics interpretation — We now show how ra- 
dius dynamics explains the Q{t) behaviors presented above. 
We will use the convenient formulation of BEC radius dy- 
namics, Eqs. (|T|) and (|2]). Fig. [3] (top row) shows the ra- 
dius of a 2D BEC evolving as a function of time for vari- 
ous ramp shapes, for reasonably large r. In the center and 
bottom rows, we show the deviation of (7{t) from the equi- 
librium radius corresponding to the instantaneous value of 
the interaction, i?o(i) = Ro\U{t)\. For a truly adiabatic 
ramp, (j{t) would follow i?o(t) exactly; therefore the devia- 
tion /(t) = (j{t) — Ro{t) is at the heart of non-adiabaticity and 
this quantity determines the heating Q. After the ramp, the 
f{t) function is purely oscillatory; Q scales as the square of 
the oscillation magnitude. The oscillations initiated at the be- 
ginning of the ramp are of magnitude 0{t~"-) for the (t/r)" 
ramp. In the case of [A] ramps (middle row), the endpoint 
kink causes an 0{t~^) oscillation, which is parametrically 
larger for a > 1 and hence dominates the final dynamics, 
leading to overall Q{t) ~ behavior. 

We first explain the ^ scaling of oscillations initiated 
at the kink. If we neglect the smaller oscillations at t < r, 
the radius a{t)KiRo{t) at the kink t — t has "correct" value 
for U = Uf, i.e. / is negligible. However the derivative 
is nonzero, a' {t)KiR'i^[t)\t=T, which scales as ^ t^-^. Thus 
we have the following "initial" conditions at < = r+ for 
subsequent evolution: /(r) = 0, /'(t+) = cqt^^. Us- 
ing f{t > r) w psin{nt + S), these initial values imply 
p ~ T^^. This explains the 0{t^^) oscillation magnitude 
of f{t > t), and hence 0{t^'^) heating, for ramps having an 
endpoint kink. 

The oscillations of Q{t) (Fig. |2] left panels) can be ex- 
plained by relaxing the approximation /(t < t) « made 
above. The small oscillations of f{t < t) guarantee that 
a'{t ~ t) oscillates around i?g(i = t) as a function of r. This 
causes final breathing mode amplitude p to oscillate around 
its 0{t^^) value as a function of t. The f{t < r) oscillation 
strength is 0{t~") (shown below); hence the Q{t) oscilla- 
tions decay as r~^", as seen in Fig.[2]left panels. 

For smooth [B] ramps (Fig.[3]bottom), the breathing-mode 
strength (~ t^") initiated at the beginning of the ramp re- 
mains unchanged; there is no kink to abruptly create larger 
oscillations. We therefore need only to explain the strength 
of oscillations at the beginning of the ramp, where s{t/T) « 
(t/r)". For this, we rewrite Eq. ([TJ as an equation for f{t). 
For simplicity, we will write this out explicitly only in the 
limit ?7i f ^ 1, and small oscillations, f{t)<^Ro{t). (The ar- 
guments can of course be modified to go beyond the large-C/ 
restriction. Small f{t) is guaranteed for large r.) We obtain 
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FIG. 3: Top row: radius dynamics a(t) for various ramp shapes, r — 
15. Center and bottom rows: deviation f{t) from the 'instantaneous' 
ground-state radius Ro{t). 2D BEC; (Ui,Ui) = (20,21). 

with u = L//(27r)^/^. The first two terms give pure oscil- 
lations (breathing mode at fixed u); the last two terms are 
corrections due to time-varying interaction. We first treat 
a > 1 ramps. The initial conditions at f = 0+ are then 
/(O) = /'(O) = 0. With u ^ u, + {5u){t/T)°', the u" cor- 
rection is dominant compared to the u'^ correction at t <C r. 
The dominant correction terms take the form ci / for a = 2, 
and cit/r^ for a = 3. The solutions of the resulting differ- 
ential equation are sums of oscillatory and algebraic terms. 
It is straightforward to verify that the boundary conditions 
/(O) = /'(O) ~ force the oscillatory part to have coeffi- 
cients scaling as ~ t^". This explains the Q ~ t--2" ijg. 
havior for integer a > 1. The a = 1 case is sUghtly differ- 
ent. The initial condition still involves cr'(O) = 0, but since 
= 1/(^+1) has finite slope at t ^ 0+, this now 
corresponds to /'(0+) = —Rq{0^) = —c^/t. This initial 
condition leads to a purely oscillatory f{t) with amplitude 
~ r^^, which explains Q{t) ^ t^^ for a = 1. 

Discussion; Open questions — Considering various 
harmonic-trapped atomic clouds, we have presented results 
on the adiabaticity question, demonstrating universal aspects 
in the first corrections to adiabatic behavior for slow ramps. 
We have shown that the slow-ramp response is determined by 
the radius-oscillation modes common to many trapped atomic 
systems, and its exact form depends on the ramp shape. We 
have explained the universal effects using a single-radius de- 
scription of the cloud. This covers a wide range of interest- 
ing systems, but it is even more remarkable that our universal 
features extend to at least some systems which cannot be de- 
scribed by a single radius. 

We have shown that a final kink in the ramp shape plays 
a drastic role in the non-adiabatic response of compressible 
trapped systems. A recently discovered effect of such kinks 
is logarithmic contributions io Q(t) ||2[|3,01- The effect we 
have found for trapped systems (kink induces larger oscilla- 
tions overwhelming initial excitation) is quite different. 



Our work opens up several new research avenues. Like 
other "universal" results, it is important to identify the lim- 
its of validity. For example, do the same Q{t) asymptotic 
features appear in trapped systems not described by a single 
radius? Paradigm examples are Bose-Hubbard systems con- 
taining superfluid-insulator "wedding-cake" structures, and 
phase-separated imbalanced Fermi gases near unitarity. We 
have shown some examples where the same behaviors appear 
[Figs. |2c,f)], but a general understanding is lacking. It is 
likely that the dynamics of one radius-like variable generally 
dominates the extreme asymptotics, recovering our results. 

Each of the systems are of intense interest in their own 
right, and understanding less universal features in parameter 
quenches is important for the individual systems, especially 
with growing experimental interest and capabilities for study- 
ing ramps and quenches. Ramps in the trapping frequency 
should also induce radius oscillations, but details might be 
different from interaction ramps. The case of anisotropic har- 
monic traps also remains an open issue. One might specu- 
late that one of the trapping frequencies dominate the extreme 
asymptotics, but the intermediate-r region might show inter- 
esting interplay of the several frequencies and associated radii. 
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